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REPORT No. 323 


FLOW AND FORCE EQUATIONS FOR A BODY REVOLVING IN A FLUID 


IN FIVE PARTS 

By A. F. Zahm 


SUMMARY 

This report, submitted to the National Advisory Committee jor Aeronautics Jor publication, is 
a slightly revised Jorm oj U. S. Navy Aerodynamical Laboratory Report Ao. 380, completed for the 
Bureau of Aeronautics in November, 1928. The diagrams and tables were prepared by Mr. F. A. 
Louden; the measurements given in Tables 9 to H were made for this paper by Mr. R. H. Smith, 
both members of the Aeronautics Sta^. 

Part I gives a general method for finding the steady-fiow velocity relative to a body in plane 
curvilinear motion, whence the pressure is found by Bernoulli^ s energy principle. Integration of 
the pressure supplies basic formulas for the zonal forces and moments on the revolving body. 

Part II, applying this steady -flow method, iinds the velocity and pressure at all points of the 
flow inside and outside an ellipsoid and some of its limiting forms, and graphs those quantities for 
the latter forms. In some useful cases experimental pressures are plotted for comparison with 
theoretical. 

Pari III finds the pressure, and thence the zonal force and moment, on hulls in plane curvi- 
linear flight. 

Part IV derives general equations for the resultant fluid forces and moments on trisymmetrical 
bodies moving through a perfect fluid, and in some cases compares the moment values with those 
found for bodies moving in air. 

Part V furnishes ready formulas for potential coefficients and inertia coefficients for an ellipsoid 
ami its limiting forms. Thence are derived tables giving numerical values of those coefficients for 

a comprehensive range of shapes. 
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FLOW AND FORCE EQUATIONS FOR A BODY REVOLVING IN A FLUID 


PART I 

INTRODUCTION 


Steady-flow Method, — In some few known cases one can compute the absolute particle 
velocity q' at any point (x, y, z) of the flow caused by the rotation of a body, say with uniform 
angular speed U, in an infinite inviscid liquid otherwise still. Thence, since q' is unsteady at 
[Xj ijj z)j the instantaneous pressure there is found by Kelvin’s formula pvlp = — d<pldt — g'V 2 , Pv 
being the supervacuo pressure there, and <p the velocity potential. 

Otherwise superposing upon said body and flow field the reverse speed — about the same 
axis, gives the same relative velocity q but which now is everywhere a steady space velocity. 
In the body’s absence the circular flow speed at the radial distance R would be qo = If 

the fixed body’s presence lowers the speed at (x, 7/, z) from qo to 3, it obviously begets there the 
superstream pressure 


P=9p(<Zo"-3-) 


( 1 ) 


or in dimensionless form, a being some fixed length in the body, 


V 

Ipa^a- 

2 


R^ 




iur 


The present text finds p by this steady-flow method only, and applies it to streams about various 
forms of the ellipsoid and its derivatives. 

The superposed circular flow, q^^ -~ 9 .R= —b\p bRj has the stream-function 

( 2 ) 

which, for rotation about the 2 axis, plots as in Figure 4 . This flow has no velocity potential, 
since b^jbRy^O. 

General Formulas for Velocity Components, — In plane flow,^ as is known, a paiTick' 
at any point (x, y) of a line s drawn in the fluid has the tangential and normal velocity com- 
ponents 

_^ip _d\f/ ,, 


1 Thia velocity entails the centrifugal pressure Po=pii*/2V2at all distances, R— from the rotation axis of the circular stream, here assumed 

to be constrained by a coaxial closed cylinder infinitely large. To the dynamic pressure pc+p may also be added any arbitrary static pressure 
such as that due to weight or other impressed force. 

* At any surface point of the body q is the velocity of wash or slip, whether the body moves or not; it is ?'i— the difference of the tangciuial 
space velocities of the fluid and surface point. If the body is fixed q"i=0, q=q‘t. 

> Plane flow, viz two-dimensional flow, literally means flow in a plane; the term applies also to space flow that is the same in all parallel islunes. 
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where 5s, 5n are elements along the line and its normal. As usual, are reckoned positive 

respectively along 5s, 5n positive; e. g. Figure 2. The components along x, y are 


^ = = 

In solid flow (3), (4) still hold for and further w^'b<pl'bz. In 
eral, + + At any point of a surface drawn in the 

rl IS t ft Wen in the olane of a and On- All these velocities are referred 



general, + At any point of a surface drawn in tne / 

fluid qt is taken in the plane of q and qn- All these velocities are referred /h^. 

to fixed space. / / \ 

Surface Velocity, — A fixed body in any stream, since = has / j 

the surface flow velocity g = which put in (1) determines the surface ^ 

pressure . . . ^ ^ 

At any surface point of an immersed moving bod}" is the same for ^ 
body and fluid, hence is known from solid kinematics. Thus, if the body “ 
is any cylinder rotating as in Figure 1 , any '‘kw cylinder having an 

^ gular speed fi about any axi.s 

g„= -nfl AElAs^^R sin (0- *3) = Q/i, =n(mx-Zy) (5) 

d/J/(l« = mx-/if, 1 , m 

where the symbols are as defined in Figures 1, 2. being daection cosines of the 

More generally, for any surface with velocities Q,,, % about the ^ 

axes X, y, z, “ * '*““*■ 9-^W^-nwdn. 

' ^ ' At any surface point Qn is the 

/ , /T \r\ y /W.V.. 7 ^,\n same for body and fluid; qt 

(/„ — (ny ~ 4 - (mx ly)^z \ ) different except at points of no 

1 * 1^10 \ slippage 

w'here Z, w, R are the direction cosines of the surface normal, as in (.loi). 

If at the same time the body has translation components, Z7, F, IF along x, y, 2, (6) must be 

increased by Z LX m F + n IF, giving 

y 

q^ = l(U-r z 9 .^-yQ,) + m{V+xn^-zQz) + n{W+yil^-xQ„) -- 0 ) 

X But (5), (6), (7) express q„ only at the model’s surface. 

\ Equations (1) to (7) obtain whether the fluid is inside or out- 

\ side the body. 

b ^\\ Zonal Forces and Moments. — For any cylinder spinning 

////^ about z, as in Figure 1 or 5, surface integration of p gives, per 

/ X, 4 unit of 2-wise length, the zonal* forces and moment, respectively, 

A ly\ ^ x = fpdy Y=J-pdx iy=y^,rdr_._--(8) 

K.ov!: ..-oco^emc data for confoca. whcrc p dy, p dx are the X y components of the elementary 
ellipses. x=a' cos 17=^ cos fi; y=b' sin surfacc force p ds, and r IS the radius vector oi (x, y). lo aerive 
yj=r sin P-, ~ tan 5=tan t,-~ tan 0 = note that p ds has Components pr d^,p dr along and across 

iY »..r“sia(.-»ffa=reos(9-ffl. r. Having no moment, p r d^ can be ignored, leaving only p 
beingeccen. dr with arm r. Thus, 2 N= fp d(H), whicn vanes as the area of 
tricity Of o6 graph of p vcrsus 

A surface oj rotation about x, spinning about its z axis, has zonal forces 

x = ffp dy dz Y=ffp dx --(9) 

•A rone is aay part'of tba surface bounded by two parallel places; in this text they are assumed normal tor, and the loae has the bounding 
planes 1 = 0 , z = ±ii; in Part III other planes are used ; e. g. i=ij. i = a. 


y 

A 


// 

y \ 

///rA 

s\\ 

/ 

\ \ 

A 

A 

-i\ / ac// a a 


Figure 2. — Geometric 'data for confocal 
ellipses. i = a' cosi 7 =rcos 0', y=b' sin 

i} = T sin 0\ ^ tan 5=tan tan 0 = 

V "i ’ (0—0)> cos (9-^0). 

/=oe = a' e\ «= Vl-6*/a* being eccen- 
tricity of ab 


201 



REPORT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


If ds,, ds^ are elements of its lines of meridian and latitude, as in Figure 3, the moment about 
z oi p ds,f dsw is r dr ds« in the plane co = 0, and p r dr dsw cos u = p r dr d 2 = dA^ for any 
meridian plane; hence the zonal moment is 


N^fPrdr. 


where P = 


J p d 2 = dF/dx, is the y-wise pressure-force per unit length x-wise.® Thus, as for 

-Zo 


(8), iV varies as the area of the graph of P versus P. Also one notes that 


Y=fPdx 


p coso? dcj (lOi) 


Since p is synunetrical about the x axis, Z = 0 = A"; viz, the assumed zone is 

not urged along y, z or about x, y, z. In general, X is not zero for such a zone, but is zero for 

the whole model. The zonal F, X are zero for stead}' 
y rotation about 2 in a frictionless liquid, because p is 

^ symmetrical about the x axis; but are not so in a viscid 

/ \ fluid, nor for accelerated spin in a perfect fluid. 

F / y \ For trisymmetrical surfaces we note also: If the 

/ \^ zones were formed by planes normal to 2 , zonal A' 

/ \ V M would be zero for motion about 2 ; zonal X in general 

/ \ P. y\ \ ^\\ zero; e. g., for a viscid fluid. Similarly for zoms 

lis N ^ with faces normal to y. 

i / ^ /a , 

\/ By (10) the bending moment about the 2 ordinate 

/ z ra 

\ y in the plane 2 / = 0 is J P r dr. This is zero for a frie- 
ze tionless liquid; for a viscid fluid it increases with length 

/z of zone. 

PiGURK 3 . ^Geometric data for prolate spheroid. i=acos Ii^ addition to the pressure forccs and moments just 
T,;y”6sin.,cos«-=r sin/3cos«; z=b sin ij sin c.>=r sin considered, duc to rotation about 2 , a viscid fluid exerts 

^ysliifeZycated Surface friction symmetrical about the 3 axis, but not 

treated here. 

For any surjace S, clearly (9) still holds and (10) can be generalized to the usual form 

X^ J'J'p(xdx~ydy)dz (IO2) 

Geometrical Formulas. — Alost of the surfaces treated in this text are members of the 
confocal ellipsoid family 


dr. This is zero for a fric- 


X 6"+ \ P 




whose semi axes are a' = + X, etc. The following known properties are needed. 

The distance from the center to the tangent plane at the point (x, y, 2 ) of a'6'c' is 




<-> 


The direction-cosines of the normal to said plane are 


j ~hf^x fh^y ~tt^z 

Ij Tfly 71 — ”72 ^ y' 2’ ^7“2- 


* The radius of the latitude circle is denoted by z#=> 
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The partial derivatives of X are 




dx dy 

More generally for any surface /(x, ?/, z) —0, one knows 


-( 14 ) 




-(13,) 


and the distance from the origin to the tangent plane at (x, y, z) is 

Ji 2 = lx + 7ny-{-nz^r cos y (12) 

7 being the angle between the radius vector r and the normal. 

Conventions. — In all the text x, y, 2 have the positive directions shown in Figure 3, 
a? also have the x, y, z com- y 



as previously stated. For a closed surface is positive outward and ds is positive 
in the direction of one walking on the outer surface wdth n on his left. 

The word “displaced fluid/^ used in treating the motion of a submerged body, usually 
means fluid that would just replace the body if the latter were removed. 



REPORT No. 323 


FLOW AND FORCE EQUATIONS FOR A BODY REVOLVING IN A FLUID 


PART U 

VELOCITY AND PRESSURE 


(A) BODIES IN SIMPLE ROTATION 

Elliptic Cylinder. — For an endless elliptic cylinder, of semiaxes a, 6, c (=«>), rotating 
about c with angular speed He in an infinite inviscid liquid, otherwise still, one knows ^ 

= ~^m' c^ca'b' sin 2t) = -~m' cos 2rj (15) 


V 



the geometric symbols being as in Figure 2. For any outer confocal a'y the potential coefficient 
has the constant value 

m',^ia + by{a'-y)/2a'b\a'-tb') (16)‘ 


On the model's surface a'~a, h' = h; {a‘^ — b^)/2ab. 

The equipotential lines on either surface ab or a'6' are its intersections with the corre- 
sponding family of h 3 ^perbolic cylinders ~ const. Normal to the equipotentials 

are the streamlines ^ = const. Graphs for ^^ = 0, 0.2, 0.4, etc., are shown in Figure 5 for a model 
having a/5 = 4. The^" are instantaneous streamlines, and form with the model a constant 
pattern in uniform rotation about c in said infinite liquid. 

At any outer confocal a'5' the velocity components are, if K^7n\a'¥9.„ 


f d<p n 

2 


. drj 

= ^ = Ksm2,5-^= 


-q', tan 2?). 


.(i: 


» Proofs of (15), (23), (29), (40) are found in books; e. g., Lamb §S 72, 106, 110, 115, 5th ed., except that Lamb reverses the sign of 

* Equivalent to (16) is wi'e-* ( ^ f > e, e‘ being the eccentricities of ab, a’b'. On ab this becomes ra',“fVVl— <*. See (49) 
\e «+Vl— e'*' 2vl— 

for the six potential coeCadents m,, mt, me, m\, m't, m'«, in the value of for more general motion. 
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\vl)ore dr 7 /ds = 1/a' — cos^?;, as one easily finds. Alternative to (17) are 

q*t= -m\Ucr cos {0 + 0) q\= -q\ tan 2r] (17i) 

Thus for t?^0, 45°, 90° (17) and (17i) give g',/Oc= — m'ca', o, m' . At the niodeTs surface, 

where (a^— Z>^)/2a6, (17i) become 

q't^ cos {0 + 0) q\=-Q,r sin {0-0) (17a) 

the latter being hfic? as in (5). 

Where q'i = 0j or cos r} = ll-\^, viz, at the stream poles, clearly x = y = 

x^“ = aV/2 (18) 

a rectangular hyperbola. (18) is the instantaneous polar streamline, e. g., Figure 5, orthogonal 
to all the confocal ellipses. Its asymptotes are +x\ its vertices are at x= :tae/V2; it cuts 

each ellip.se where xjy^a'iy, viz, on the diagonals of the circumscribed rectangle. For an 


endless thin plate of wddth 2a the poles are ht y = 0j x^ ±a/V2. 

Superposing on the body and fluid, and using (2), changes (15) to 

\t^ = ^{r^ — m' ca'y cos 2r?)flc (1^1) 

Its graph, with A^ = 0.2, gives the streamlines in Figure 6 for the flow 11^ = "1 round a fixed 


cylinder having a/b = 4. About the point (0, 1.45) in Figure 6, is a whirl separated from the 
outer flow by the streamline ^ = 4.25. This line abuts on the model at the inflow points f, i; 
spreads round it and emerges at the outflow points o, o.^ The streamlines for an endless thin 
rc(“tangle having 5 ==0, e =lj are similar to those of Figure 6, but infinitely crowded at the edges. 
The superposed particle velocity — contributes to (17i) 

q'\= cos {0-0) = sin {0-0) = (20) 

also q" tan {0—0). Adding (17i) and (20) gives the components qt = q\ + q'\t + 

of the resultant flow’ velocity at any field point. One notes that (20) are the reverse of qi, 
in Figure 1. 

Ill particular qn^O on the fixed model and x, y axes; hence there 

(//al2c= cos (0 + 0) + cos {0—0)] qjqo = m\ cos {O + 0)+cos {0—0) (21) 

Tims qlqo^ 1 + on the x axis; 1 — ?n'c on the y axis; and 1 at oo where w'^ = 0. The dashed 
line in Figure 6 gives qja^c^ — {\ — m' c)yja for points on the y axis; it crosses y at the whirl 
ceiiliT where ^|[ = 0, viz, w’here m'c=l. By (16) w'c^l for the surface of any model having 
(i;h>\ + ^[2; and there is no whirl if a/6<Cl + V2. Figure 7 shows for the surface of a 

model having a/6 = 4, m'c= (a^ — 6^)/2a6 = 15/8. 

Putting of (21) in (li), where rVa^ = cos^i;/cos^/3, gives 

p/^pa^l2^c = (1 “ cos {O + 0) + cos {O — 0)Y) cos^ 77 /cos- d ---(22) 

which is graphed in Figure 7 for a model having alb=4, 

integrating p/~pa^Q‘c, as in (8), gives for an inviscid liquid T=0 = A"; X?^0. Figure 7 
(ielineates X for this case. 

Ttie points t o are identical with those in Figure 5; viz, where the slip speed g in (21) is zero; they are called stop points, stagnation points, etc. 
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For the surface of an endless flat plate (6 = 0, c= «) fixed in the stream — clearly 
m'c = a/26 and generally r cos (0 — j8) =0; hence (21) gives 

T 

5/aQc= — cos (0 + ^)— - sin 9 cos tj cot 2r? (2li) 

which equals — oo^ o, 1/2 for = 45^^, The flow resembles that in Figure 6; it has 

twin whirls abreast its middle, stop points at x= ia/V2, and infinite velocity at the edges. 
Putting in (li) r^x and qo= —xQc gives the platens surface pressure 

= (1-cot^ 2v) cos“ r, (22,) 




IiGURB 6. — Streamlines about endless elliptic cylinder fixed m infinite inviscid liquid rotating about its long axis with uniform angular speed — 
shows ^ — y fl (r*— m'i a* b' cos 217) with increments A = n= — 1. Dotted line portrays z-wise speed on v axis 

which equals —1/4, 1/2, — <» forx = 0, ±a/V2, ±a; viz, for t; = 90°, 45°, 0, etc. 

Prolate Spheroid. — For a prolate spheroid, of semiaxes a, 6, c, rotating about c with 
speed ilc in an infinite inviscid liquid, 


— m' cQtXy^ sin 2 v cos i 
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the geometric symbols being as in Figure 3. For any outer confocal spheroid a'b'c' (23) has 
the known constant potential coefficient 


3 , 1+e' e'^ 

3 2 N 1 1 ~r € n t 

2e(2-*)logi-^-6 + i_ 


ee 


(24) 


€ being the eccentricities of oh, a h . Table gives surface values of 7 ti\ for various shapes 
of prolate spheroid. 

In the yzy zx planes (p = 0; in the plane, where cos a; = 1 


sin 2t) \P= cos 2t? 


(23,) 


which, except for have the same values as (15), entailing the same polar streamlines (18). 
The equipotentials on a'5'c' are its intersections with the family x?/= — const. 

At any point (x, y, z) on a'6'c' the orthogonal velocity components are by (23) 


f d ip 

^ dn 


, b(p d 77 


y doj 

ds„ 


(■25) 


&n, «s„ 5s„ denoting line elements along the normal, meridian, and circle of latitude, as in 
Figure 3. Since g'„ is absent from (1), we shall not need it; we merely note that on the model's 
surface it is sin (d—fi) cos u. By geometry drj/ds, = r cos {8 + p)la'b' cos 2ij,* da)/ds„ = 
1/6' sin 7?; hence 

q\= —m'cQcT' cos (^-h/S) cos oj (/u = cos /3 sin o> (25, ) 

For w = 0, o' differs only by m' ^ from (170 for an elliptic cylinder; also r cos (3 = x 

q^u> = m\ xQc sin w = 0, m'c^Qc for O7 = 0, tt/2. 

Superposing — on the above system adds to (25,), as easily appears 

2 n ~ sin (^— ^) cos cj cos (6 — ^) cos w 2 "„ = cos ^ sin w (21)) 

At the now fixed surface and on the x, y axes Q_n^0 = q' n + hence summing (250, (26) 
gives there 

Qv~ “[w^cCoe (^-f/3)-f cos (0-/3)] cos cos w| 

2«= (1 +m'c)^^r cos /3 sin co=g^sin cj j' ' 

Thus for 07 = 0 clearly qlqo = m'c cos (^ + |3)-f cos (<9-/3), differing from (21) only by for 

o) = x/2, q/qQ = ~ (1 +m '0, a formul a like that for a negative flow qo across a cylinder; for = 

90°, 45°, q = qm (S^-f-Fc.)- On the x axis 2 / 2 o=l + m',; on the y axis q/qo^l~ m\>0 

everywhere, hence no whirl centers on y. 

Figure 8 shows Iq/aUci on the meridians a; = 0, ± 45°, ± 90° of a fixed spheroid witli a;b - 4. 
Distributions symmetrical with these occur on the opposite half of the surface. Notewortliv 
is g for w = ± 90°. By (27) it is g = ± (1 4- m'<,)f2cX; hence the straight-line graph in Figure s. 

Figure 8 shows also, for these meridians, the pressure computed wuth the working foriniilii, 
derived from (li), (27). 

V r> 

y- —A cos“ CO + ^ sin^ o) \' 2 s) 


* E. g., by (23) 2 ^ ^ = cos 2r, . which gives ^ in (25). Also directly g'^ 


-rn'eSlef COS cos U. 
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where cos (0 + /3)+cos id-^)Y) cosVcos^/?, - m',(2 + m',) cos^r;. Here m% = 

.689 by Table IV. The crosses and circles, giving experimental values taken from Keference 3 
show good agreement with (28) for a considerable part of the surface. For cos co = 0, ^oc5oc.r; ' 
or the graph is parabolic. 

Integrating^, as in (9), (10), gives for an in viscid liquid Y^0-=N,X^0. Figure 8 portrays 
X computed from theory and experiment. 

Ellipsoid, — For an ellipsoid, of semiaxes a, 6, : along x, ?/, rotating about c with speed 
ilc in an infinite inviscid liquid, otherwise still, 


<p= 

which for any outer confocal ellipsoid a'5'c', has the constant potential coefficient 


(29) 


m'c^ C{p-a) 


a^-¥ 

2{a^-¥)- (d^ + ¥) (/3o ~ «o) 


(30) 


the Greek letters being as in Part V. Surface values of are listed in Table IV. 

By (29) the equipotential lines on a'6'c' are its intersections with the hyperbolic cylinder 
family xy= ~ <plm\nc = const. The orthogonals to <p const, at the surface a'&'c' are the stream- 
lines there. These by (31) are parallel to x where x^O; parallel to y where y = Q; normal to 
^ where z = 0. The same obviously holds for spheroids and other ellipsoidal forms. 

In the xy plane the flow has the polar streamlines (18); also it has there 


(p = --^m ' sin 2?; 


\p= - cos 2t) 


( 29 :) 


whence the streamlines in that plane are plotted. The form of (29i) is like those of (15) and 
(23i), for the elliptic cylinder and prolate spheroid, entailing similar expressions for the velocity 
and pressure in the plane-flow field z = 0. 

For the general flow the velocity components at a'5'c' are by (29) 

«'= + «'= + w' = - (31) 

and those due to the superposed velocity — QcR = 5'oj are 

w"=-0 (32) 

whence the resultant velocity and pressure may be derived for all points of the flow field about 
the ellipsoid fixed in the steady stream —UcR. In forming the x, y, s derivatives of 7n\ one 
may use the relations (14) and (72). 

Everywhere in the planes x = 0, y = 0j the resultant velocities are respectively, by (31) 
and (32), 

q = u= {l~m\)Ucy q==v = - (1 + (33) 

while in the plane z = 0, q can be found as indicated for an elliptic cylinder. (33) apply also 
to the elliptic cylinder and prolate spheroid previously treated, and to all other forms of the 
ellipsoid fixed in the flow ~ 

(B) BODIES IN COMBINED TRANSLATION AND ROTATION 

Most General Motion. — The most general motion of any body through a fluid may 
have the components 77, F, W along, and Oa, about, three axes, say a, 6, c. The entailed 
resultant particle velocity q' at any flow point is found by compounding there the individual 
velocities severally due to 77, F, F, Oj, and computable for an ellipsoid by formulas in 
Reference 2 and the foregoing text. 

Yawing Flight. — In airship study the flow velocity q' caused by a prolate spheroid in 
steady circular flight is specially interesting. Let the spheroid's center describe about 0,. 


210 



FLOW AND FORCE EQUATIONS FOR A BODY REVOLVING IN A FLUID 


9, a circle of radius na, with path speed nail. Then if a is the constant yaw angle of 
Mttack, the component centroid velocities along a, ft, and the steady angular speed about c are, 
respectively, 

U = nail cos a V— nail sin a = (34) 

If, now, velocities the reverse of (34) are imposed on the body and fluid, ?n = 0, and the 
siii face velocity q on the fixed spheroid has in longitude and latitude the respective components 

g, = (1 +A:a) sin 0 - {l+ 1 ct)V cos 6 cos co- [m'^ cos (0 + /3) + cos { 6 -^)]ilcr cos w] 

= (1 +Z:j) F sin W+ (1 + m\)ilcr cos jS sin w j 

where positive flows along ds are, respectively, in the directions of increasing tj, w, as in Fig- 
ure 3. The terms in V, are known formulas for translational flow, e. g., Reference 2; the others 
are from (27). Hence q^ then p is found for any point (/3, w) on the spheroid.^ If fic is negli- 
gible, q = q sin €, where q^= (l+JCayU^+ {I and e is the angle between the local and 
|)olar normals, as proved in Reference 2. 

Figure 02 portrays, for specified conditions, theoretical values of p/^pQ^j Q being the path 
sp(‘od spheroid's center; it also portrays pl^pQ^ for the model in rectilinear 

motion, with Q= U. The difference of p/^pQ^ for straight and curved paths, though material, 

is loss than experiment gives, as shown by 9a. Fuller treatment and data are given in 

K(‘forcncc 3, 

The forces X, Y and moment N, for any zone, may be computed as before; but for the 
whole model they are more readily found by the method of Part IV. Zonal Y and N for a 
hull form are found in Part III. 

The first of (35) applies also to an elliptic cylinder, with cos m\= {a~ — ¥)j 2 ah. 

Fixed in a flow — Uy — F, —ilc, it has the surface velocity 

q= {\ + h/a) V sin e ~{l+a/i)V cos d- + <^) + c®® (»“ /3) (36) 

For an endless flat plate 6^0, cos B = hja. sin 0 cot t?; and the last term of (36) may be rewritten 
by (21 1); thus (36) becomes 

g = ( U— V cot 1) — ailc cos tj cot 2 r}) sin 6 (37) 

These two values of q with (li) give the pressure distribution over an elliptic cylinder or flat 
plate revolving about an axis parallel to its length or fixed in a fluid rotating about that axis. 

Thus an endless plate of width 2a, revolving with angular speed 12, path radius na, and 
inc idence a, as in Figure lOi, has by (37) the relative surface velocity, viz, slip velocity 


qfail={n cos a — n sin a cot 17 — cos tj cot 2 r}) sin 6 (38) 

and since sin 1, qo^=U^-\- (F-f acl2)^ = a^‘^(n^ + 2ri sin a cos 77 Fcos^tj), (1) gives 

2>/4pa^l2^ = ’l^ + 2n sin a cos t/Fcos ^r) — n^ (cos a — sin a cot cos 77 cot 277)^ (39) 

j n 


For ?x = 3, a = 30®, Figure IO2 delineates the distribution of slip velocity q/ail on both sides of 

the plate; IO3 that of the pressure p!^ pa’^il^ on its two faces. This pressure integrated over 

the plate’s double surface gives F=0, as may be shown. The dashed line in Figure IO3 is the 
pressure-difference graph whose integral for 77 = 0 to tt is also zero. The resultant forces X, Y 
and moment N for such a plate are found in Part IV by a method simpler than surface inte- 
gration of the pressure. 

* Here again q is the slip speed of the flow at any point of the body’s surface, and depends only on the relative motion of body and fluid. 
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Flow Inside Ellipsoid. — At any point inside an ellipsoid with speeds Z7, V, W, ^2*. 
Qcy along and about a, h, filled with inviscid liquid otherwise still, 


ij ^2 — q2 ^2 — q2 q2 

V? — Lx -f Vy 4- 2 + _|_ ^2 LlaV^ ^ -I- 6^ - 


■ (40) 



Figure 9. — Prolate spheroid in steady yawing flight. (1) Defines 
velocity conditions; (2) delineates theoretical pressure dis- 
tribution; (3) experimental pressure distribution for Q=40 feet 
per second. In (2) and (3)» full lines indicate rectilinear, 
dashed lines curvilinear motion 


Figure 9 (continued).— For conditions (1), (4) delineates pressure load (Jtr unit 
length; (5) the zonal force; (6) the zonal moment. In (4) the full and dotted lines 
give theoretical values from equations (oi), {bi); the dashed line, experimental values 
from reference 3. (5) is obtained by planimetring (4); (6) by planimetring (j) 


whose coefficients are constant for the whole interior. Hence the components of the particle 
velocity g are 


bx 


— a^ — 


.(41) 
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and like values for v, w found by permuting the symbols. If the fluid were solidified any 
particle would have 

etc., etc. (42) 

Thus when an ellipsoid full of in viscid still fluid is given any pure translation its content moves 
as a solid; but when given pure rotation each particle moves with less speed than if the fluid 
were solidifled, since the fractions in (41) are less than unity. 

For velocities ?7, F, of the ellipsoid 



Figure 10.— ‘Endless flat plate revolving about axis parallel to its length, in infinite Inviscid fluid. (1) Defines conditions; (2) delineates 
relative velocity qla Q of fluid; (3) pressure a and pressure differenoe ^ p/-^ p 0^ on two faces of plate 


for which w = d(p/bz = 0. For this plane flow (4) with (43) gives 

Vy- (44) 

whence the streamlines may be plotted. In particular if the model has simple rotation Uc, 

+ 

= -2^2^2^/i2,=-const (45) 

and the interior streamlines are hyperbolas, as in Figure 5. 

Adding (2) to ^ in (45) gives the steady flow 

^ = + - (46) 

hence the streamlines lie on the elliptic cylinders 

4* = (a^ -h h^)f^c^ — const (47) 

104397—30 28 
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By (46) 5 = + which put in (1) gives at (x, y), since ^o= 

« -r, A(«y+_*y)_ ^^8, 

where 2?n = P5oV2. Here is the centrifugal pressure due to the fluid’s peripheral velocity 
qoy and p is the pressure change due to go “5, 2 being the relative velocity of fluid and container. 
In a like balloon hull g would quickly damp out, leaving only as the dynamic pressure. At 
the ends of a, h, c, respectively, (48) gives 

Pn-p_ 46^ 4a^ 

Pn {a^ + hY 


For large a/h the first is negligible, the second approaches 4, giving ~3pn = — 1.5pQ/6^ as 
the temporary dynamic pressure drop inside the hull at the end of 6. Experimental proof 
would be interesting. 

Potential Coefficients. — An ellipsoid of semiaxes a, h, c along x, ?/, z, when moving 
through an infinite inviscid liquid, otherwise still, with velocities U, F, IF, ^b, Uc along and 
about the instantaneous lines of a, 6, c, begets the known velocity potential 

(p= —maUx — mbVy~mcWz—m\^ay2-7n'bilb^-m\n^y ( 49 ) 

the six potential coefficients m being constant over any outer confocal ellipsoid a'6'c'. Their 
values for abc are given in Tables III, IV. Alternatively (49) can be written for this surface 


<P=-~'ka Ux - lc„Vy - Icc ^ 


the Ics being the more familiar inertia coefficients defined and tabulated in Part V. Of the six 
potential coefficients in (50) the first three are the same as the inertia coefficients A’ a, kb. 
the last three are greater except when cjb or ajc or ft/a is zero. Thus, if ft/a = 0 the last term of 
(50) is — Jc'cU^yy which is the potential on the outer surface of an elliptic cylinder (a^ oo ) rotating 
about c. Everywhere inside of it the potential is ficXy, as (40) shows. 

For the flow (40) textbooks give the inertia coefficients 



which are the squares of the potential coefficients. One notes too that the ratios of like terms 
in (40), (50) equal the ratios of like potential coefficients and like inertia coefficients, which 
latter in turn are known to equal the ratios of like kinetic energies of the whole outer and inner 
fluids, if the inner moves as a solid. 

Relative Velocity and Kinetic Pressure. — When a body moves steadily through a 
perfect fluid, otherwise still, the absolute flow velocity it begets at any point (x, y, z), being 
unsteady, is not a measure of the pressure change there. The relative velocity is such a measure. 
To find it we superposed on the moving body and its flow field an equal counter velocit}", thus 
reducing the body to rest and making the flow about it steady. The same result would follow 
from geometrically adding to said absolute flow’ velocity the reversed velocity of (x, y, z) assumed 
fixed to the body. In particular this process gives for any point of the body’s surface the wash 
velocity, or slip speed, w’hich w’ith Bernoulli’s principle determines the entailed change of surface 
pressure. Conversely, if the pressure change at a point is known or measured, it determines 
the relative velocity there. In hydrodynamic books the above reversal is used commonly 
enough for bodies in translation. In this text it is employed as well for rotation; also for 
combined translation and rotation. However general its steady motion, the body is steadily 
accompanied by a flow’ pattern wffiose every point, fixed relatively to the body, has constant 
relative velocity and constant magnitude of instantaneous absolute velocity and pressure. 
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FLOW AND FORCE EQUATIONS FOR A BODY REVOLVING IN A FLUID 


PART m 

ZONAL FORCES ON HULL FORMS » 

Pressure Loading. — For a prolate spheroid abc with speeds Z7, F, Figure 9i, or fixed 
in a stream — f/, — F, — Qc? (35) gives at {x, ?/, z) on abc the relative velocity 

(f — qj A — B cos 03 + C cos^ co 

A, By C being constant for any latitude circle. In forming this equation one finds 
= 2(1 + /(:«) U sin B{ (1 + F cos d + [m\ cos (0 + )3) + cos (0- 
etc., for Aj C. In the body’s absence said stream has, at said point (x, y, 2 ;), 

"" ( “ bJ -\-y^cY + { — F“ x^cY=Ai — Bi cos oj -f (7i cos^ co, 
where co alone varies on the latitude circle. Its radius being yo = ZQj makes y — yo cos co, 

Bl=2UZo^lcy 

etc., for Aif Ci. Putting g, go in (1) gives the surface pressure 

pl.5p = qQ^~q^= (^1 — ^4) + {B~Bi) cos co+ (f7i — C) cos^ co. 

r2ir r2ir 

By (lOi) the loading per unit length of x is, since I cos co = 0= I cos^ co, 

F/.5p= — f p cos codco = — (F— Fi)zo r cos^ codco= ^ 7 t(F — Fi)2q (a) 

•^pjo Jo 


Ay Axj Cy Ci vanishing on integration of p. Thus, finally, 

PI.5pQ^=-^{B^Bi)zolQ^ (a,) 


P having the direction of the cross-hull component of g? at co = 0. 

One notes that g«(oc sin^ co) contributes nothing to B or the integral in (a) ; viz, the loading 
P is unaffected by q^, and depends solely on g„ the meridian component of the wash velocity. 
Also for ^ = 0 and tt, B — Bi = 0 = P. 

In Figure 94 the full line depicts (aj) for the spheroid shown in 9i, circling steadily at 40 
feet per second. The theoretical dots closely agreeing with it are from Jones, Reference 3, as is 
also the experimental graph. Beside them is a second theoretical graph plotted from Doctor 
Munk’s approximate formula derived in Reference 8 and given in the next paragraph. But 
that Professor Jones omitted some minor terms in his value of p^ his theoretical Pj.5pQ* should 
exactly equal (ai). His formula, derived by use of Kelvin’s pv/p = ^ — gV2, can best be studied 
in the detailed treatment of Reference 3. 

In Reference 8 Professor Ames derives Munk’s airship hull formula 


P 

r^pQ- 


sin 2a 


d5^ 2 d 
dx ^Fdx 


(xS)y 


1 This part was added after Par^ I, II, IV, V were typed; hence the special numbering of the eauations. 
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S being the area of a cross-section; R the radius of the path of the ship’s center. This was 
assumed valid for a quite longish solid of revolution; for a short one it was hypothetically 
changed to 

={Jc ,- fc J sin 2a ^ {xS) (b) 

Applying this to a prolate spheroid we derive the working formula 


P 


-Lx-Mx^^-N. 


where the constants for a fixed angle of attack are^ 


• 7T sin 2a, 


¥ 

i/=3Fc 

a" 


2x 

cos a, 


(b.) 


gCOS a. 


Plotting (bi) for the conditions in 9i gives the dotted curve in 94 . It shows large value.< 
of PI.5pQ^ for the ends of the spheroid, where (aO gives zero. To that extent it fails, though 
with little consequent error in the zonal force and moment at the hull extremities. It has th< 
merit of being convenient and applicable to any round hull whose equation may be unknown 
or difficult to use. 

Zonal Force.—Aii end segment of the prolate spheroid, say beyond the section x = 
bears the resultant cross pressure 




which with the resisting shear at Xi must balance the cross-hull acceleration force on the seg- 
ment in yawing flight. For the whole model (bO with (c) gives F-0, which is not strietb 
true for curvilinear motion; but (ai) with (c) gives the correct theoretical value of Y, am! 
agrees with (67). 

In Figure 9s graphs of Yj.bpQ\ for the values (aO and (bO of P, are shown beside tho-o 
derived from Jones’ experimental pressure curve. Since F is proportional to the area of 
segment of the graph of P, it can be found by planimetering the segment or by integrating Pd.r. 

Zonal Moment.— The loading P exerts on any end segment, say of length a-x, tlie 
moment about its base diameter z 

A'. = JVdx 


which can be found by planimetering the graph of Y. Figure 9j delineates so derived from 
the three graphs of Y. They show the moment on the right hand segment varying in length 
from 0 to 2a; also on the left segment of length from 0 to 2a. The resisting moment of the 
cross section must balance Nz and the acceleration moment of the segment. 

Correction Factors. No attempt is here made to deduce theoretically a correctifci 
factor to reconcile the computed and measured p. In Reference 3 Jones shows that the tlu - 
retical and e.xperimental graphs of P/.5pQ- have, for any given latitude *,> a/2, the san e 
difference of ordinate whatever the incidence 0</a</20°. Thus the ordinate difference fouml 
for the zero-incidence graphs, when applied to the theoretical graph for any fixed 0</a<^2i> , 
determines the experimental one with good accuracy. Such established agreement in loadiim' 
favorably affects, in turn, the graphs of F, the transverse force and moment on any emi 
segment of the spheroid. 


‘ From the meridian curve 


=*1 


o_ 

dx a* yo ’ 


, L dS „ iitjQ 0* 

iryo*; hence -=2jryo - -2 jt x, which put in (b) leads to (bi). 
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FLOW AND FORCE EQUATIONS FOR A BODY REVOLVING IN A FLUID 

PART IV 

RESULTANT FORCE AND MOMENT 

Body in Free Space. — Let a homogeneous ellipsoid of semiaxes a, ft, c move freely with 
component velocities Vy Wy respectively along and about instantaneous fixed space 

ixes Xy yy z coinciding at the instant with a, 6, c. Then the linear and angular momenta referred 
oXy yy z are 

m]U TYiiV rriiW Axp C]T (52) 

being the body’s mass, A^y B^y C\ its moments of inertia about a, 5, c. If, now, forces Xiy 
)\y Zi and moments Zi, J/i, are applied to the body along and about x, y, 2 , they cause in 
the vectors (52) the well-known change rates 

mi{u~rv+qw)^ Xi Aip— {Bx— Ci)qr = Li 

mi{v — pw + ru)^Yi Biq— {Ci — Ai)rp= Mi (53) 

7rii{w — qu + pv) = Zi Cir— {Ai — Bi)pq= Ni 

which apply to any homogeneous solid symmetrical about the planes ahy he, ca. 

For motion in the ah plane; viz, for w, p, q = 0; (53) give 

Xi = 7rii{u~rv) Yi = mi{vi-ru) Ni = CiV (54) 

and for uniform revolution about an axis parallel to 2 , as in Figure 11, viz, for tt, v, f = 0, (54) 
become 

Xi=—mirv Yi = miru Ni=^0 (55) 

where now Xi, Yi are merely components of the centripetal force mir^Ju^ i-v^, whose slope is 
Yi/Xi = —ujv. Also if (2^ + is the path velocity of the body’s centroid, h its path radius, 

r ^ Qjh is the angular velocity of h and of vector tyiiQ. 

Reactions of Fluid. — If external forces impel the ellipsoid from rest in a quiescent fric- 
tionless infinite liquid, with said velocities u, y, w, p, q, r, they beget in the fluid the corresponding 
linear and angular momenta 

Icamu lct,mv hemw k^aAp k\Bq k\Cr (56) 

w here m is the mass of the displaced fluid, and A, B, C its moments of inertia about a, h, c. 

One calls kc7n the apparent additional masses’^; ¥ aA, ¥ i,By ¥ cC the “apparent 

additional moments of inertia,” of the body for its axial directions; because the fluid’s resistance 
to its linear and angular acceleration gives the appearance of such added inertia in the body. 
The six ¥s are called “inertia coefficients,” and are shape constants. Values of them are 
given in Tables III, VI, VIII for various simple quadrics. 

The component flow momenta (56), like (52), are vectors along the instantaneous directions 
of a, h, c; viz, along x, y, z\ hence their time rates of change must equal the forces and moments 
w hich the body exerts on the fluid; viz, 

X==m(JcaU-kirv-\-lccqw) L = ¥ aAp~ (Jc\B-¥ cC)qr~ (ki, — ]cc)mvw 

Y^m{kt,v~kcpw + ]caru) M=¥ i,Bg— {¥ cC--¥ aA)rp— (Jcc — ka)mivu (57) 

Z = m{lccW — lcaqu-\-kt,pv) N = ¥ cCr— (¥ aA — ¥ tB)pq~ {ka—kb)muv 

‘ These Dew meaoiDgs of u, v, w, jT, q, r are assigned for convention's sake and for convenience. 
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all written from (53) on replacing its momenta by those of (56), and adding vector-shift terms. 
Thus the vector 'k^rnw shifts with speed v entailing the change rate Iccmw.v of angular momentum 
about Xy while Ici^mv shifts with speed w entailing the opposite rate — Ic^mv.w. Their sum is 
{kc — lc^)rmw. Permuting these gives for the ?/, 2 axes (k c)mwfiy — When the 

h^s are equal the vector-shift terms vanish, as for said free body, or for a sphere, cube, etc., in a 
fluid. The fluid reactions are (57) reversed. (57) ai)ply also to fluid inside the trisymmetrical 
surface. 

If the angle of attack is a^tan~*v/u, we may write in (53), (57) 

r^Qjh u==Q cos a v^Q sin a uv = ^Q^ sin 2a (58) 



Figuke 1 1 .— Momenta and forces for free body in uniform circular motion. Centripetal force, /i?i = 77iQr=* 
TnQVA, has slope — u/r, r being angular speed about 0 


Of special aeronautic interest are (57) for plane motion, such as in yawing airship flight. 


for which Wy p, 5 =0, giving 

X = m{kaU~ki,rv) Y=7n(kt,v + karu) N^=k' cCr+ Qc^ — ka)7nuv (5t0 

Thus for uniform circular flight 

= ~k i;mrv Y = k a w ru N = {ku~ka) muv ( 6 i ) ) 

which are the analogues of (55) for the free body. Or in notation (58) 

A'= sin „ cos « ‘V = sin 2a (til) 

T being the volume of the model. 
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As shown in Figure 12 (60) give the resultant force and slope 

= mr^llca^u^-^^bV^ YIX=—^ cot a= — cot ^ (62) 

Kj, 

also R and N at the origin arc equivalent to a parallel force R through the path center 0, along 
a line (called the central axis of tlie force system) whose arm and intercepts are 

N/R = h sin {0~ a) x=^lsecp y = Z cosec 0 (63) 



Kkilke i2.— Momenta and forces for symmetrical body in uniform circular motion through frictionless inflnite 
liquid otherwise at rest. Whole hydrodynamic force, R = mry/ka^ has slope ~ka ujkh v. Yaw 

moment N—ikh—ka) inuv = {kb — ka) t ^ sin 2 a, t being volume 


For steady motion (60) show that the body sustains no force in pure translation (r-0); 
no force nor moment in pure rotation (w, ^; = 0); no moment in revolution about a point on 
X or ?/; viz, for u = 0, or = For given u, v the moment is the same for revolution as for pure 
tiiiuslation. The forces result from combined translation and rotation; the moment from 
translation oblique to the axes a, 6, irrespective of rotational speed. 

Combination of Applied houcEs. find the whole applied force constraining a bodv 
to uniform circular motion in a perfect fluid (55), (60) may be added, or graphs like those of 
Figures 11, 12, may be superposed. For an airship having = (55), (60) give 

X~ — (\-\-lCi>)mvr y ~ {].+lca)mur N= (Jci,~lca)^uv (64) 

* Writing rj=»rQ.m cos^a-j-fti^ sin^a we may call it the centripetal force of the apparent mass m sin"^ for the body direc- 

tion I'vf Q. 
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where Z = + etc. Figure 13, compounded of Figures 11, 12, shows that a submerged 

plane-force model, revolving uniformly about its path center, may have as sole constraint a 
single force R through that center, and outside itself ; that is attached to an extension of the 
model. Such conditions appear commonly in vector diagrams of aircraft. The line of S, so 
defined, is the central axis of the force system. 

Hydrokinetically Symmetric Forms. — Equations (56), (57), for trisymmetrical shapes, 
apply also to others having hydrokinetic symmetry. Examples of these are: All surfaces of 
revolution, axially symmetric surfaces wEose cross sections are regular polygons; torpedo forms 
svmmetrically finned, etc. All these figures, as has been known many decades,^ have three 



Figure 13.— Composition of forces on symmetrical borly in uniform circular motion through frictionless infinite liquid 
otherwise at rest. Resultant of centripetal and hydrodynamic forces, R = m r p’, has slope 

’ Figure 13 is 11 and 12 compounded 

1+jCfc V 


orthogonal axes with origin at the body’s impulse center,^ such that if the body, resting in 
quiet sea of perfect fluid, is impelled along or about either axis it begets in the fluid a line.u 
or angular momentum expressible by a vector along that axis. 

Examples. — We may apply (60) to some simple cases interesting to the aeronaut!* d 
engineer. 

(1) For an endless elliptic cylinder in uniform yawing flight, as in Figure 12, m^Tpah p- r 
unit length, and by comparison with Table VIII lca = hja, kf, = alh] hence by (60) 

* See Reference 7. 

* I. e.» the point of intersection of k^m U, k^mV, kcm W; it may be found as in the last paragraph of Part V. 
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N= Tr{a^ — ¥)p.uv = 7r(a^ — 6^) 2 sin 2a. 


( 65 ) 


A'' = — ira^prv Y= ir^pru 

The resultant force Tcpr^jd^v^ + has the slope — &%/a^z;= — 5^/a^.cot a; the central axis is 
through the path center; X is the same as for a round cylinder of radius a; Y the same as for 
one of radius h. For a good elliptic aircraft strut a/6 = 3; hence A’/F= — 9?;/^^ = —9 tan a; 

A=87r&^pU2? = 87rJ^. 2 -sin 2a, By (65) N is the same for all confocal elliptic cylinders, since 


ii^—¥ is so. 

* 

If a — h^ as for a round strut, A^=0, /? = Tra^prQ^ and coincides with the body’s previously 
found centripetal force to which it bears the ratio m/mj. 

If 6 = 0, as for a flat plate, (65) become 


X = ~ ira^prv F = 0 


N= TTO^puv = ira^ 


pQ^ 

2 


sin 2 a 


(66)* 


The equivalent resultant force ira^prVf with slope F/A=— 0, runs through the path center 
parallel to x. If r = 0, the plate has pure translation, with forces A, F=0, and moment 
X=Tra^puVf a well known result. X in (66), being the same as in 
65), is independent of the strut thickness 6. 

(2) For a prolate spheroid, of semiaxes a, 6, 6, in uniform 
yawing flight, m = 4/3.7rpa6^, and Icaj are as given in Table III. 

Thus for alh = 4, Jca, ^* = 0.082, 0.860; hence by (60) 

X= ~3.6a6Vy Y=0'.M34a¥pru N=3.26a¥puv.{67) 

(3) For an elliptic disk of semiaxes a, 6, c, moving as in Fig- 
ure 14, Table VIII gives = ^rpabyE; hence by (57) the forces 
and moment are 


4 a 


Y= ~!ccmpw= ~ ^-^.TTph^.pw 

L = Jc pn.VW = ~^,7rpb^.VW 


Z = 0 

-( 68 ) 



Figure 14. — Thin elliptic wing moving par- 
allel to its plane of symmetry through a 
perfect fluid 


the other pertinent terms in (57) vanishing, as appears on numerical substitution. Here 
/t = £'(fl,|), sin^0= (a^-¥)/a^; also i = sin 2a. Compare (68) with (66), calling b the 

width in both. 

Theory Versus Experiment. — In favorable cases the moment formulas of Part IV 
i*(‘cord fairly well with experiment, as the following instances show. For lack of available data 
Mic force formulas for curvilinear motion are not compared with experiment. 

(1) By (65) an endless elliptic strut with a=l/3 foot, 6 = 1/12 foot, c = 5 feet, held at a 
licgrees incidence in a uniform stream of standard air at 40 miles an hour, for which pQ^l2^ 
4 093 pounds per square foot, sustains the yawing moment per foot length 


V=7r(a'-6^) 


pQ^ 

2 


sin 2a= 1.3392 sin 2a lb. ft. 


-(69) 


This compares with the values found in the Navy 8 by 8 foot tunnel, as shown in Table IX 
i;dred from Figure 15. The agreement is approximate for small angles of attack. The model 
w iis of varnished mahogany, and during test was held with its long axis c level across stream, 
and with two closely adjacent sheet metal end plates, 2 feet square, to give the effect of plane 
tl'kW. 


Equations (66) were published in Reference 5 as the result of a special research to determine the fluid forces and moment on a revolving plate 
Ii: I he present text they follow as corollaries from more general formulas. 
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(2) By (66) an endless thin flat plate of width 2a = 5/12 feet, similarly held in the same 
air stream, has per unit length the moment 


sin 2a = 0,5581 sin 2a lb. ft_ 


This is compared in Table X and Figure 16 with the values found in the Navy 8 by 8 foot 
tunnel. The flat plate was of polished sheet aluminum 3/32 inch thick, wuth half round edges 
front and rear. 

Again for an endless flat steel plate 5.95 inches wide by 0.178 inch thick at the center, with 
its front face flat and back face V-tapered to sharp edges, Fage and Johansen, Reference 6, 




Figure 15.— Theoretical and experimental moment Figure 16.— Theorectical and experimental moment about long axis of endless 

about long axis of en<lless elliptic cylinder. Width rectangular plate. Width 5 inches, air speed 40 miles per hour. Correction 

S inches, thickness 2 inches, air speed 40 miles per factor «=0.860 

hour. Correction factor k = 0.912 

found, at 50 feeT per second and 5.85^ angle of attack, iV= 0.125 pound foot as the moment 
per foot run about the long axis, computed from the measured pressure over the median section. 
By (66), a thin flat plate would have 

iV=7ra- • ■ sin 2a = 0.1931 X2.9725X0. 2028 = 0. 116 lb. ft. 

whicli is 7 per cent less than 0.125 found with their slightly cambered plate. 

(3) An elliptic disk 3/32 inch thick with a, fe = 15, 2.5 inches, w^hen held as a wing in the 
Nav}" 40-mile-an-hour stream, had the moment L versus angle of attack a shown in Figure 17 
and Table XL For this case 

siir6l= (a2~62)/a' = 875/900, 0 = 80^-24', £:= 1.03758. 

Also in (68) a = 5/4 feet, & = 5/24 feet, = 4.093; hence 

^ ^ 0.8963 sin 2a lb. ft (71) 
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which gives the theoretical values 
at small incidences. The disk 
edges without any rounding or sharpening. 

(4) For a wooden prolate spheroid 
24 inches long by 6 inches thick, 
carried as in Figure 12 round a circle 
of radius /i = 27.96 feet to the models 
center, Jones, Reference 3, found at 
40 feet per second the values of N 
listed in Table XII. For this case 
Table III gives A:* — = 0.778, and 
(61) gives 

pQ^ 

2 


in Figure 17 and Table XL The agreement is fair 
as tested was of sheet aluminum cut square at the 


\-.30 


K % Theorefica! 


(Atj — ^a)r- • sin 2a = 0.388 sin 2a. 





1 1 [— — r 

5 " 14 '' 

Angie of a ffack^ cC 


18" 


These values appear from Table XII 
not to accord closely with the experi- 
mental ones. 

Correction Factors. — Figures 
15, 16, 17 portray experimental - 5 ° 
moments, at small angles, as accu- 
rately equal to the theoretical times 
an empirical correction factor k. 

Thus amended (61) gives for the 
experimental moment 

• sin 2 a. 

For the given elliptical cylinder 
>c = 0.912 with “8°<a<6°; for 
the endless plate x = 0.860 with 
— 6 °<Ca<C 6 °; for the elliptic disk 
0.887 with — 5°<a<4°. In such 
cases one should expect to find the 

actual air pressure nearly equal to ^ V-.20 

the theoretical over the modeFs figure 17.— Theoretical and experimental moment about long axis of elliptic disk, 
forward part, but so deficient along ^ 

the rear upper surface as to cause a 

defect of resultant moment. No effort is made here to estimate it theoretically, nor to de- 
termine it empirically for a wide range of conditions. 

The measurements shown in Table X, for the flat plate, were repeated at 50 and 60 miles 
an hour without perceptible scale effect. 
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PART V 


POTENTIAL COEFFICIENTS, INERTIA COEFFICIENTS 


Green’s Integrals. — The foregoing text employs Green’s well-known integrals, which 
for the ellipsoid abc may be symbolized thus: 




.(72) 


where a' = -y^+ 6'= + > are semiaxes of the confocal ellipsoid a'Vc\ The integrals 

have the following values, Reference 4: 

a = AQ>^-c^)[F{e, <p)-E{d, ^)] 




-(73)^ 


w^here 


A = - 


2abc 


sin^^ = 




(T — r 

X" 


-(74) 


-(75) 


'{a^~¥W-¥)^o}-e 

and the elliptic integrals are 

-^(^j ¥^) ^ y'(I “ sin^0 sinV)~^^dv’ E{Qj Si\~ sin^^ sinV)^dv>- 

Numerical values of F{B, ^), E(B, ^), a, 0, 7 are given in Tables I, II for X = 0 and various 
ratios a/6, 6/c; viz, for various shapes of the ellipsoid a6c. For — tt/2 one writes F{B^ tp) — 
E(By ip) = Ey by convention. 

Potential Coefficients. — For motion (49) the ellipsoid ahc has the potential coefficients 
known from textbooks. 

m . 

” 2-fio 


, <?( t -/ 3 ) . ^ b^-c^ 

, H(a-y) ^ 

o = ^Tzrr~^Z:\ where 


mc = 


2-70 


2H— (Ofl — 7o) 

/(/3-a) 

■.27-(0 o-«o) 


= ^here 7=“^^ 

C ^ WIICIC 1 (X^^¥ 


-(76) 


rWa, Wft, Tfic being for translation along a, 6 , c and m'&, m% for rotation about them, and 
«o, /3o, To being (73) for X = 0; viz, for a', 6 ', c' =a, 6 , c. Surface values of (76), viz, for 7 = 
«oj 00 ; 7 o are given in Tables III, IV. For fluid inside the ellipsoid the potential coefficients are 
as in (40) and given numerically in Table V. 

Inertia Coefficients. — From (76) are derived the conventional linear and angular inertia 
coefficients 

^ay 'kty Icc^rria, TYic Jc\, ¥ tj ¥c = Gm\, Hm\y Iwf c (77) 

for the ellipsoid moving through or containing liquid, as in (40), (49). Surface values are 
given in Tables III, VI, VII. 

1 (73) satisfy the known relation a-k-^-\-‘t = 2 abcja‘b'c\ as appears on adding. 
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Limiting Conditions. — In sofne limiting cases, as for c = 0, or a==5, etc., (73) may become 
indeterminate and require evaluation, as in Reference 4. In such cases the formulas in Table 
VIII may be used. For c==0, entailing zero mass and infinite ICa Tc\y lc\y one may use in (57) 
the values of Iccm, ¥ aA^ ¥ given at the bottom of Table VIIL 

Physical Meaning of the Coefficients. — The tabulated potential coefficients, put in 
(40) or (49), serve to find the numerical value of the potential <p^ or impulse ~pv> per unit area, 
at any point (x, y, z) of an ellipsoid surface.^ Integration of over any surface, as explained 
for j> in Part I, gives the component linear and angular zonal impulses. So, too, integration of 
— Pv?^„/2, where is the normal surface velocity at (x, y, 2 ), gives the kinetic energy imparted 
to the fluid; and integration of the impulsive pressure — pb<p!'bt gives the impulsive zonal 
forces and moments. One finds pd^/d^ for (40), (49) by using with them the specified density 
p, accelerations Z7, V\ IF, (lay fic, and tabulated potential coefficients for the given semiaxes 
df 5, c. 

Thus putting ~p<pc, ~p<p'c for p in (9), (lOi), and integrating over the whole ellipsoid 
surface, easily gives the fluid’s linear and angular momenta 

ICc.mW ¥c.CUc (78) 

where m W, C^c are respectively the linear and angular momenta of the displaced fluid moving 
as a solid with velocities Tf, The like surface integration of - p^c?n/2 gives, as is well known, 

¥c^C^\j2 (79) 

where m TP/2, are the kinetic energies of the displaced fluid so moving. Each inertia 

coefficient therefore is a ratio of the body’s apparent inertia, due to the field fluid, to the like 
inertia of the displaced fluid moving as a solid. 

By (49) the potential coefficients due to velocities W, are 

Ulff— (pci TVs TYl c~ ^ cj^^y 

The first is the ratio of the outer and inner surface potentials due to TV at any point z on the 
ellipsoid dbc] the second is the ratio of the potentials due to 12^ at (x, y), respectively on the 
outer surface of that ellipsoid and inside the cylinder of semiaxes 00 , 5, c. 

One notes that the momenta (78) times half the velocities give (79); also that the time 
derivatives of (78) are the force and moment Z, N =lCcmiV^ ¥ cCtlcy as in (57) for the simple 
2 -wise motions, TV, (l. 

For any axial surface, say of torpedo form, moving as in Figure 12, the ratio —¥ cC^J^bT^V 
is the distance from the arbitrary origin 0i to the impulse center O 2 , or center of virtual mass. 
This may be taken as origin, and if the body’s center of mass also is there Figures 11, 12 can 
still be superposed as in Figure 13. In the same way are related the acceleration force and 
moment ¥cCUcj thus illustrating the doctrine that the motion of a hydrokinetically sym- 

metric form in a boundless perfect fluid, without circulation, obeys the ordinary dynamic 
equations for a rigid body. 


Aerodynamical Laboratory, 

Bureau of Construction and Repair, U. S. Navy, 

Washington, D. C,, December 17 j 1928. 

* This impulse is imparted hy the moving surface to the fluid, otherwise still; the fluid in turn tends to impart to the body the impulse 
per unit area at (i, p, z). 
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CHIEF SYMBOLS USED IN THE TEXT 


a, b, C- - 
a'j b'f c' 

€, e' 

n; hi, h 2 
I, m, n. _ 

S, Sif, Suj _ 
X, y, 
r, oj_ ^ 

77, e 


M, V, 10 

Qi, Qn 

q„ q 

U, V, w 

u, V, ir 

p, q, r 

Sla, ^b, 

^ 

m«, mb, me 

?n'a, m' by vi'e 

0=^lu^v^'+^v^'.. 

AuBiyCi 

A, B,C 

miy m 

P, r 

p, pn 

Xi, r., Z,; Ri 

X,Y,Z; R 

Li, Ml, Ni 

L, M,N 

hay hb, kc 


GEOMETRICAL 

Semi axes of ellipsoid abc. 

Semiaxes of confocal ellipsoid a'b'c'. 

Eccentricities of ellipse ab and its confocal a'b'; ae — a'e' = ^la^—b'^ 
Normal to ellipse ah', distances from origin to normal and tangent. 
Direction cosines of normal n to any surface. 

Length along any line; lengths along meridian and circle of latitude. 
Cartesian coordinates; also coordinate axes. 

Polar coordinates of prolate spheroid abc. 

Eccentric angle of ab, inclination to x of normal to ab. 


KINEMATICAL 


Component velocities of fluid parallel to x, y, z axes. 
Component velocities of fluid parallel to tangent and normal. 
Resultant velocity of fluid before and after disturbance. 
Component translation velocities of abc parallel to a, b, c 
Component translation velocities of abc parallel to a, b, c 

Component rotation velocities of abc about a, b, c 

Component rotation velocities of abc about a, b, c 

Velocity potential, stream function. 

Potential coefficients for abc with velocities u, v, w or U, V, W. 
Potential coefficients for abc with velocities p, q, r or i}a, 1^6, 
Resultant velocity of abc. 

DYNAMICAL 


Alternative symbols. 


Moments of inertia of rigid body about its axes a, b, c. 
Moments of inertia of displaced fluid moving as a solid. 

Mass of body, mass of displaced fluid. 

Density of fluid, volume of model or displaced fluid. 

Pressure of fluid moving, pressure on coming to rest. 
Component forces applied to free rigid body; resultant force. 
Component forces exerted by body on fluid; resultant force. 
Component moments about a, b, c applied to rigid body. 
Component moments about a, 6, c exerted by body on fluid. 
Inertia coefficients for abc moving parallel to a, b, c in fluid. 
Inertia coefficients for abc rotating about a, b, c in fluid. 
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TABLE I 

ELLIPTIC INTEGRALS F{6, <p), E{9, 


[Defined in eq. (75), Part V) 


a/C 






bfc 






1 


3 

4 

5 

6 


8 

9 

10 








Fie, ip) 





1 

0.00000 










1 

2 

1.31696 

1.04720 









i 

3 

1.76305 

1. 43870 

1.23095 









4 

2.06412 

1. 71374 

1. 48399 

1.31814 








5 

2.29319 

1.92798 

1.68471 

1.50687 

1. 36940 







f) 

2.47903 

Z 10413 

1. 85188 

1.66560 

1.52053 

1.40332 





7 

2.63508 

2.25400 

1.99520 

1. 80281 

1. 65204 

1. 52959 

1. 42745 



• 

s 

2.77024 

Z 38432 

Z 12075 

1. 92379 

1.76856 

1.64194 

1. 53595 

1. 44550 




y 

Z 89035 

Z 49971 

Z 23255 

2. 03191 

1. 87318 

1. 74.321 

1.63405 

1.54065 

1. 45948 



10 

2.99638 

Z 60288 

Z 33303 

Z 12855 

1.96804 

1. 83534 

1. 72357 

1. 62768 

1.54419 

1. 47063 



•0 


CO 

CO 

CO 

CD 

1 00 

i 

«» 

00 









E{e, ^) 





1 1 

0.00000 



1 







! 2 

.86603 

1. 04720 









3 

.94277 

1. 07024 

1.23095 

! 








4 

.96822 

1.06091 

1. 18103 

1 1.31814 








5 

. 97975 

1.05019 

1. 14337 

1. 25126 

1. 36940 







6 

.98597 

1.04146 

1. 11604 

! 1. 20294 

1.29996 

^ 1. 40332 


1 




. 98972 

1.03472 

1.09589 

1. 16833 

1. 24893 

i 1.33574 

. 1. 42745 





8 

.99214 

1.02946 

1.08071 

! 1. 14185 

1. 21035 

! 1.28451 

I. 36317 

1. 44550 




1 9 

.99378 

1. 02529 

1.06894 

' 1. 12136 

1. 18040 

1 1.24464 

1.31304 

1.38483 

1. 45948 



10 

.99496 

1. 02195 

1.05960 

1. 10516 

1 1. 15669 

1 1. 21297 

1. 27310 

1 1. 33642 

1. 40240 

1. 47063 


CO 

1.00000 

1.00000 

1.00000 

1 1.00000 

1.00000 

1.00000 

1. 00000 

1.00000 

1.00000 

1.00000 












, . 



— .. 


I The integrals in this table are culled from L. Potin's Formules et Tables Numerique. 

TABLE II 

GREENES INTEGRALS ao, /3o, 70 


[Defined in eq. (73), Part V] 


Q/C 






bjc 





j 

1 2 

3 

4 

5 

6 

7 

8 

9 

10 

00 1 


ao 

1 

0.66667 

i 









[ 

2 

.34713 

0. 47280 1 









1 

3 

.21751 

. 31265 i 

0, 36460 








: 

4 

. 15092 

.22474 

.26820 

0. 29636 







! 

5 

. 11171 

.17064 

. 20719 

.23189 

0. 24951 






1 

6 

.086527 

. 13471 

. 16584 

. 18769 

.20336 

0. 21541 






7 

.069266 

. 10950 

.13629 

.15541 

. 16970 

. 18079 

0. 18950 




, 

g 

.056894 

.091037 

. 11435 

. 13135 

. 14426 

. 15440 

. 16254 

0. 16914 



1 

9 

. (M7710 


, 097571 

. 11276 

. 12448 

. 13378 

. 14132 

. 14757 

0. 15271 



10 

.040637 

0 

.066203 

0 

.084381 

. 097957 

0 

. 10872 

0 

. 11728 

. 12428 
0 

. 13010 
0 

. 13500 

0 

0. 13920 
0 

0 







Bo 





i 

1 

0.66667 











2 

,82643 

0.47280 










3 

.89127 

.534Z3 

0.36460 









4 

. 92459 

.56964 

.39662 

0.29636 








5 

.94418 

, 59182 

.4I8M 

. 31587 

0. 24951 







(1 

.95678 

.60693 

. 43307 

.32965 

.26265 

0.21541 






7 

.96538 

.61775 

.44413 

.34083 

.27275 

. 22477 

0. 18950 





g 

.97154 

. 62577 

.45260 

.34912 

.28071 

.23234 

.19654 

0. 16914 




9 

97619 

.63184 

.45913 

. 35569 

.28709 

.23847 

.20237 

.17458 

0. 15271 



10 

97972 

. 63659 

.46437 

,36109 

.29233 

.24354 

.20725 

, 17927 

. 15712 

0. 13920 


CO 

1.00000 

.66667 

.50000 

.40000 

.33333 

.28572 

.25000 

.22222 

. 20000 

. 18182 

0 







yi 






1 

0.66667 











i 2 

.82643 

1. 05440 










1 ^ 

. 89127 

1. 15312 

1. 27078 









4 

.92459 

. 1. 20572 

1. 33518 

i 1. 40726 








5 

. 94418 

1. 23752 

1.37478 

1. 45223 

1.50098 







6 

. 95678 

1. 25835 

1.40110 

1. 48267 

1, 53401 

1. 56918 





i 

7 

.96538 

i 1. 27273 

1.41956 

1.50377 

1, 55754 

1. 59443 

1.62100 





8 

.97154 

1 1.28320 

1. 43306 

1. 51953 

1. 57504 

1. 61325 

1.64091 

! 1. 66172 



1 

' 9 

. 97619 

1 1. 29109 

1. 44329 

1.53154 

1. 58844 

1. 62775 

1. 65630 

1. 67784 

1. 69457 



10 

. 97972 

1. 29720 

1, 45125 

i 1.54094 

1. 59895 

1. 63917 

1.66846 

1. 69062 

1. 70787 

1. 72160 


' CO 

1! 00000 

1 1. 33333 

1.50000 

^ 1.60000 

1.66667 

1. 71429 

1.75000 

1. 77778 ' 

1.80000 

1. 81818 

2. 00000 

' 



. 

... 




- - - 



- ■ 

— 
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TABLE III 

POTENTIAL COEFFICIENTS m», m.* FOR ELLIPSOIDS IN TRANSLATION 

(For outer surface of a b e) 

[Defined in eq. (76)] 



b/c 

0/c 

' 

2 

3 

4 

5 

6 

1 

■ 7 

8 9 f 10 ® 



77»o 

1 

0.5000 







1 1 

2 

.2100 

0 . 3096 






1 ! 

3 

.1220 

.1853 

0 . 2229 






4 

. 08162 

, 1266 

. 1549 

0. 1740 





5 

.05916 

.09328 

. 1156 

. 1312 

0. 1425 




6 

.04522 

.07222 

.09042 

. 1036 

. U 32 

0 . 1207 



7 

.03588 

.05792 

. 07313 

.08425 

.09272 

.09938 

0 . 1047 


8 

.02928 

.04769 

.06064 

.07029 

.07774 

.08366 

. 08846 

0.09238 

9 

.02444 

.04008 

.05129 

. 05975 

.06637 

. 07169 

. 07603 

. 07966 0 . 08267 

10 

.02074 

.03423 

.04405 

.05150 

. 05748 

.06229 

.00626 

.06958 . 07239 0.07481 


0 

- 

' 0 

0 

0 

0 

0 

0 

0 0 0 0 



1 

0.5000 







1 

2 

.7042 

0. 3096 







3 

.8039 

.3645 

0 . 2229 






4 

.8598 

.3983 

.2474 

0. 1740 





5 

.8943 

.4203 

.2643 

. 1876 

0.1425 




6 

.9171 

,4357 

.2764 

. 1974 

.1512 

0. 1207 



7 

.9331 

.4469 

.2855 

.2054 

.1579 

. 1266 

0. 1047 

i 

8 

.9447 

.4554 

.2925 

.2115 

. 1633 

. 1314 

. 1090 

0.09238 

9 

.9535 

.4618 

.2980 

.2163 

. 1678 

. 1354 

. 1126 

. 09564 0 . 08267 i 

10 

. 9603 

,4669 

.3024 

.2203 

.1712 

.1387 

. 1156 

.09846 . 08526 I 0.07481 


1.0000 

.5000 

.3333 

.2500 

.2000 

.1667 

.1429 

.12500 .11111 .10000 0 







m . 



1 

0.5000 







1 

2 

.7042 

1.115 







3 

.8039 

1 . 362 

1.743 






4 

.8598 

1.518 

2 . 008 

2 . 374 





5 

.8943 

1 . 623 

2 . 199 

2. 651 

3.008 




6 

.9171 

1. 697 

2. 339 

2.866 

3.292 

3. 642 


i 

7 

.9331 

1.750 

2 . 446 

3. 030 

3.520 

3.931 

4 . 277 

1 

8 

.9447 

1 . 790 

2 . 528 

3 . 163 

3.706 

4 . 171 

4 . 570 

4.912 

9 

.9535 

1. 821 

2.593 

3.269 

3.860 

4 . 373 

4 . 819 

5 . 208 5 . 548 

10 

. 9603 

I . 846 

2 . 645 

3.357 

3.987 

4 . 543 

5.032 

5.465 5 .^ 6 184 

» 

LOOOO 

2.000 

3.000 

4.000 

5.000 

6.000 

7.000 

8.000 9.000 10.666 oa 1 


— 

— 


....... 

— .-. 

, 



1 


• These have the same values as the inertia coefficients fc#, kb, k,. 
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TABLE IV 

POTENTIAL COEFFICIENTS m'a, m'c FOR ELLIPSOIDS IN ROTATION 

(For outer surface of a 6 c) 

(Defined in eq . ( 76 )] 


o/c 

1 


2 

3 

4 

1 ^ 

6 /c 

' ^ 


8 : 

1 “ 

10 











1 

0 


i 

! 









2 

0 


0 . 5643 










3 

0 


.6390 

1 . 045 









4 

0 


. 6768 

1 . 135 

1.499 








5 

0 


.6989 

1 . 190 

1.596 

1.943 







G 

0 


.7125 

1. 225 ' 

1.663 

2 . 042 

2.380 






7 

0 


.7211 

1.249 

1 . 7 a 5 

2 . 113 

2 . 481 

2.813 





8 

0 


. 7270 

1.266 

1.738 

2. 165 

2.556 

2.915 

3 . 245 




9 

0 


.7315 

1 . 278 

1.762 

2 . 20.5 

2.615 

2.995 

3 . 348 

3 . 675 



10 

0 


i .7348 

1.288 ; 

1.780 

2 . 236 

2.660 

3.058 

3 . 430 

3 . 778 

4. 103 



0 


! .7500 

1. 333 1 

1 

1.875 

i 2.400 

! 1 

2.917 

3 . 429 

3.937 

4 . 444 

4 . 950 



m'i 


1 

i 0 










2 

i - 0. 3990 

- 0 . 5643 










3 

5819 

-.8853 

- 1.045 









4 

1 -.6888 

- 1 . 104 

- 1 . 349 

- 1.499 








5 

-. 7581 

- 1.264 

- 1 . 588 

- 1.800 

- 1.943 







6 

1 -.8058 

- 1.384 

- 1 . 780 

- 2.052 

- 2.243 

- 2.380 







1 -.8402 

- 1. 476 

- 1.935 

- 2.264 

- 2.504 

- 2.680 

- 2.813 





8 

8659 

- 1 . 548 

- 2 . 062 

- 2.445 

- 2 . 732 

- 2 . 948 

- 3 . 114 

- 3 . 245 




9 

-, 88.57 

- 1.607 

- 2 . 168 

1 - 2.600 

- 2.931 

- 3. 188 

- 3.388 

- 3.547 

- 3, 675 j 

- 4 . 103 


10 

-.9013 

- 1.654 

- 2 . 257 

1 - 2.734 

- 3. 107 

- 3 . 402 

- 3 , 637 

- 3 . 825 

- 3. 978 ! 



- 1.0000 

- 2.000 

- 3.000 

I - 4.000 

- 5.000 

- 6.000 

- 7,000 

- 8.000 

- 9.000 

- 10.000 i 



m'a 


7 " 








i 


0.3990 

0 









.5819 

0.1556 

0 








.6888 

.2420 

0.08332 

0 







.7581 

.2969 

. 1359 

0 . 05193 

0 




: 

.8058 

.3350 

.1719 

.08705 

0. 03549 

0 




.8402 

.3627 

. 1981 

,1134 

.06127 

0. 02569 

0 



.8659 

. 38.36 

.2181 

. 1330 

.08081 

.04529 

0 . 01951 

0 


.8857 

.3998 

.2336 

i . 1484 

. .09610 

.06058 

.03486 

0.01527 

0 

.9013 

.4127 

.2460 

! . 1608 

. 1084 

.07291 

.04721 

.02770 

0. 01238 0 


j 1.0000 

.5000 

.3333 

i ,2500 

.2000 

. 16667 

. 14286 

. 12500 

.11111 0.10000 

0 







. , 



— 



— — 


104397 — 30^—29 
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TABLE V 

POTENTIAL COEFFICIENTS m'a, rn't, m'e FOR ELLIPSOIDS IN ROTATION 
(For all points inside of a 5 c) 
f Defined in eq. (40)] 


d/c 


afc 

1 

2 

3 

4 

5 6 1 

7 

8 

9 

10 ® 

0 to <» 

0 

j 0, 60000 

0. 80000 

1 1 

0.88235 

0. 92308 0. 94595 

0. 96000 

0.96923 

0. 97561 

1 0. 98020 1. 00000 


m'i = // 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0 

-0.60000 


-.88235 
92308 
94595 
-.96000 
96923 
-. 97561 
-.98020 
- 1.00000 


•Same for all values of 5/c 


m'e = I 


1 

2 

0 

0,60000 i 

0 


1 







3 

.80000 1 

0. 38462 

0 

i 







4 

.88235 

.60000 

0.28000 

^ i 







5 

. 92308 

. 72414 

. 47059 

0. 21951 1 

0 






6 

. 94595 

.80000 

.60000 

.38462 : 

0. 18033 

0 





7 

.96000 

.84906 

.68966 

. 50769 i 

. 32432 

0.15294 

0 




8 

.96923 

.88235 

. 75342 

. 60000 

.43820 

.28000 

0. 13274 

0 



9 

. 97561 

.90588 

.80000 

. 67010 

.52830 

.38462 

. 24615 

0.11724 

0 


10 

.98020 

.92308 

.83486 

. 72414 

.60000 

. 47059 

. 34228 

. 21951 : 

0. 10497 

0 

<» 

1.00000 

1.00000 

1.00000 

1.00000 

1 

1.00000 

1. 00000 

i 1.00000 ^ 

I. 00000 

1.00000 

1.00000 
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TABLE VI 

INERTIA COEFFICIENTS 1 k'ty k'c FOR ELLIPSOIDS IN ROTATION 
(For outer surface of a 6 c) 

(Defined in eq. (77)] 


6/c 


ajc 

1 

2 

3 

4 

1 5 

! 6 

7 

8 

9 

10 

00 






A 41 “ GfJTi Q 





1 

0 




1 

1 





! 


2 

0 

0. 3386 





1 



i 


3 

0 

.3834 

0. 8359 









4 

0 

.4061 i 

.9081 

i 1. 323 








5 

0 ! 

.4104 1 

.9519 

i 1.408 

1.793 







6 

0 ' 

. 427.5 1 

.9803 

1.468 

» 1. 885 

! 2. 251 






7 

0 

. 4326 . 

. 9995 

1. 505 

1.950 

' 2. 347 

2. 701 





8 1 

0 

.4362 ' 

1.013 

1. 533 

1.999 

j 2.418 

2. 799 

3. 145 




^ I 

0 

. 4389 1 

1.023 

1. 555 

2. 035 

2. 473 

2. 875 

3.245 

3. 585 



10 1 

0 

. 4409 ! 

1.030 

1. 571 

2.064 

1 2. 516 

2.935 

3. 324 

3.686 

4.022 


00 

0 i 

. 4500 i 

1.067 

1.654 

2. 215 

2. 759 

3. 291 

3. 816 

4 . 336 

4.852 





1 

0 









1 2 

0. 2394 

0. 3386 








3 

.4655 

.7082 

0. 8359 







4 

. 6078 

.9745 

1. 191 

1. 323 






1 5 

. 6998 

1.167 

1. 466 

1.662 

1.793 





6 

. 7622 ■ 

1.309 

1. 683 

1.941 

2, 122 

2.251 




7 

' .8066 

1.417 

1. 857 

, 2. 174 

2.403 

2. 573 

2. 701 



8 

. 8393 

1. 501 

1.999 

2. 370 

2,648 

2.857 

3. 019 

3. 145 


9 

! .8641 

1. .567 

! 2. 115 

' 2. 536 

2.860 

3. 110 

3. 305 

3. 460 1 

3.585 

10 

! .8834 

1. 622 

i 2. 213 

1 2. 679 

3.045 

3.335 

3. 565 

3. 749 1 

3.900 

(D 

1.0000 

2.000 

! 3.000 

1 

4.000 

i 

6.000 

6.000 

7. 000 

8.000 

9.000 




1 

0 











2 

0. 2394 

0 










3 

. 4655 

0. 05985 

0 









4 

.6078 

. 1452 

0. 02333 

0 








5 

1 . 6998 

.2150 

.06393 

0. 01140 

0 







6 

1 . 7622 

.2680 

. ia3i 

.03348 

0.00640 

0 






7 

1 .8066 

. 3079 1 

. 1367 

. 05758 

. 01987 

0. 00393 

0 





8 

. 8393 

.3385 ! 

. 1643 

. 07982 

.03541 

.01268 

0.00259 

0 




9 

! .8641 , 

. 3622 ' 

. 1869 

.09941 

. 05077 

.02330 

.00858 

0. 00179 

0 



10 

. 8834 

.3810 i 

.2054 

. 1164 

.06503 

.03431 

.01616 

.00608 

0. 00130 

0 

0 1 

CO 

1.0000 

.5000 

.3333 

.2500 

.20000 i 

. 16667 

. 14286 

. 12500 

. mil 

0. 10000 



1 For translation k„, kb, kt are given in Table III. 
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TABLE VII 

INERTIA COEFFICIENTS h\, k\, fc'. FOR ELLIPSOIDS IN ROTATION. 

(Inner aurface of ab c) 

(Defined in eq. (77)] 








b/c 






afc 

1 

2 


4 

5 

6 

7 


9 

10 

i ® 













0 to ® 

0 

0. 36000 

0.64000 

0. 77854 

0.85208 

0. 89482 

0. 92160 

0. 93941 

0. 95181 

0,96079 

1.00000 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

0 

0. 36000 
.64000 
. 77854 
.85208 
.89482 
.92160 
.93941 
.95181 
.96079 
1.00000 

Same for all values of bic. 




















2 

3 

4 

5 

6 

7 

8 
9 

10 

0 

0.36000 
.64000 
.77854 
.85208 
.89482 
. 92160 
.93941 
.95181 
.96079 
1.00000 

0 

0. 14793 
.36000 
.52438 
.64000 
. 72090 
.77854 
.82062 
.85208 

1. 00000 

0 

0. 07840 
. 22145 
.36000 
.47563 
.56764 
.64000 
.69699 
1.00000 

0 

0. 0481S 
. 14793 
.25775 
.36000 
.44903 
. 52438 
1.00000 

0 

0. 03252 
. 10518 
.19202 
. 27910 
.36000 
1.00000 

0 

0. 02339 
.07840 
. 14793 
.22145 
1.00000 

0 

0,01762 
.06059 
. 11716 
1.00000 

0 

0. 01375 
.04818 
1.00000 ^ 

0 

0. 01102 
1.00000 : 

0 

1.00000 1 

1 









1 

1 
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TABLE VIII 

INERTIA VALUES FOR LIMITING FORMS OF ELLIPSOIDS a>b>c 


INERTIA COEFFICIENTS FOR TRANSLATION AND ROTATION 

a/b 

Shape 

k. 

kb 

kt 

k\ 

k\ 


c=0 

1 

Circular disk 

0 

0 




0 1 

1-f- 


0 

0 


00 

00 

0 ! 


Long rectangle 

0 

0 

00 

00 

oo 

0 1 











6>c>0 




1 

J-i- 

Oblate spheroid c*/a’ 

c ce—a si 

Q-Jf 

sin~ie 

a Of — c sin-*« 
c.ce— a sin“‘e 

C'(70 

-m 

J-C»)(7o-^o')i 

0 

Elliptical cylinder 

0 

cfb 

bfc 

i (6?-c»)2 

26c 6»+c* 

bjc 

c/6 




c-6 





1 



H 

H 

0 

0 

0 

1 + 

Prolate spheroid <*=1— cVa* 

Round cylinder 


log« 

l+e 2e 
l-e l-ei 

0 




0 

, « l-2<* 

1 1 1 

0 

1 1 1 










APPARENT MASSES AND MOMENTS OF INERTIA WHEN c=0 



a/b 

Shape 

kaTn 

kbin 

ktm 

k'aA 

k'bB 

k',C 

1 

1 + 

Circular disk 

Elliptical disk 

0 

0 

0 

0 

0 

0 

> 

UfttbVE 

16 

45 

4irp ab*(a^-bi) 

oa» 

4irp o36>(a>-6*) 

0 

0- 

0 


15 {gai-b^)E-b^K 
iirp6<* 

15 (o>- 26 i)E+ 6 »Ar 

irpbi* 







* Per unit length of model. 
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TABLE IX 

LIFT, DRAG, AND MOMENT ON ENDLESS ELLIPTIC CYLINDER 

(Width 8 inches, thickness 2 inches, air speed 40 miles per hour] 


Lift 


Drag 


i Moment about long 


Angle I 
of at- 
tack u* 



' foot run 

Theoreti- 

degrees 

Pound per foot 
run 

Experi- 

mental 

cal 

N= 1.3392 
sin 2a 

-8 

-2.30 

0. 160 

-0. 335 

-0. 3691 

-6 

-1.94 

. 139 

-.254 

2784 

-4 

-1.42 

.122 

-. 170 

1864 1 

-3 

-1.11 

, 116 

-. 127 

-. 1400 

-2 

-.76 

. Ill 

-.084 

-.0934 

-1 

-.40 

. 108 

-.042 

-. 0467 

0 

0 

.106 

0 

0 

+1 

+.41 

. 108 

+.044 

+. 0467 

2 

.80 

. Ill 

.085 

.0934 

3 

1.13 

. 116 

. 129 

. 1400 

4 

1.44 

.123 

. 171 

. 1864 

6 

1.90 

. 140 

.249 

.2784 1 

+8 

+2. 16 

.165 

, +.325 

+.3691 


•As the test angles a were in part fractional, all measurements in Table IX are faired from the original graphs of lift, drag, and moment versus 
a, in fig. 15. 


TABLE X 

LIFT, DRAG, AND MOMENT ON ENDLESS THIN FLAT PLATE 

[W'idth 5 inches, air speed 40 miles per hour] 


Angle 
of at- 
tack a, 
degrees 

Lift Drag 

Pound per foot run 

Moment about long 
axis pound foot 
per foot run 

' Theoreti- 
Eiperi- cal 

mental A'‘-0.5581 
sin 2a 

-8 

-1. 345 

0. 190 

-0. 107 

-0. 1538 

-6 

-.980 

.112 

-.097 

-.1160 

-5 

-. 827 

.0816 

-.083 

-.0984 

-4 

-.614 

.0596 

-.064 

-. 0777 

-3 

-.471 

.0464 

-.060 

-.0583 

-2 

-.315 

.0360 

-.032 

-. 0389 

-1 

-. 157 

.0324 

-.016 

-.0195 

0 

0 

.0312 

0 

0 

+1 

+.155 

.0328 

+.017 

+.0195 

2 

.311 

.0360 

.033 

.0389 

3 

.471 

.0472 

,050 

.0583 

4 

.639 

.0648 

.066 

.0777 

5 

.831 

.0900 

.085 

.0984 

6 

1.018 

. 124 

.098 

. 1160 

8 

1.346 

.206 

.107 

. 1538 

10 

1. 538 

.291 

.084 

+. 1909 

12 

1.594 

.360 

.074 


14 

1.682 

.422 

.059 


16 

1.681 

.480 

.055 


+ 18 . 

, +1.530 

.542 

+, 046 
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TABLE XI 

LIFT, DRAG, AND MOMENT ON THIN ELLIPTIC WING 

[Length 30 Inches, width 5 inches, air speed 40 miles per hour] 



Lift 

Drag 

Moment about long 


axis, pound foot 

Angle 
of at- 









tack a 




Theoreti- 

degrees 

Pounds 

Experi- 

mental 

cal 

L-0.8963 





sin 2a 

-8 

-2.415 

0.426 

-0. 173 

-0. 2471 

-6 

-1.835 

,266 

164 

1863 

-4 

-1. 186 

. 169 

-.109 

-. 1247 ; 

-3 

-.886 

. 138 

-.082 

-. 0937 , 

-2 

-.667 

.116 

-.063 

-.0625 

-1 

-.294 

.106 

-.031 

-.0313 ! 

0 

-f. 005 

.103 

0 

0 


.3C5 

.106 

-f.030 

-1-. 031S 

I 

2 

.590 

.118 

.066 

.0625 

3 

.890 

.136 

.084 

.0937 

4 

1. 195 

.168 

.in 

.1247 

6 

1.861 

.265 

.166 

.1863 

8 

2.474 

.422 

.166 

.2471 

10 

1 2.865 

.667 

.134 

-f.3066 

12 

2.958 

.696 

.109 


14 

2. 892 

.798 

.094 


16 

2. 859 

.897 

.086 


18 

2.769 

.974 

.097 


-1-20 

-1-2.7^ 

1.065 

+. 095 



TABLE XII 

MOMENT ON PROLATE SPHEROID i 

[ r.ength 24 inches, diameter 6 inches, through-air speed 40 feet per second] 


Moment about minor axis, pound foot | 


Angle ;^n'5fbS-^ Found by pressure 


of at- 
tack a, ' 
degrees 


-20 

-10 


on bal- 
ance 


Recti- 

linear 


integration 


Recti- 

linear 


Curvi- 

linear 


motion motion motion 


-0. 179 
-.106 


-0. 207 

-.122 


-0. 157 
-.078 


Theoreti- 

cal 

N=0.388 
sin 2a 


-0. 249 
-. 133 


0 

— . 
0 

0 

+.021 

0 

-flO 

+.106 1 

+.122 

.127 

+.133 j 

i +20 i 

+. 179 

+.207 ■ 

+. 177 

+.249 


I Data taken from Reference 3. 
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